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ABSTRACT 

We report the first observation of Lya pumped molecular hydrogen emission lines in planetary 
nebulae. The H2 emission observed in the ultraviolet spectra of NGC 6853 and NGC 3132 can be 
explained by excitation of vibrationally hot H2 by Lya photons. Constraints are placed on the nebular 
Lya emission profile, as well as the molecular hydrogen temperature, column density and turbulent 
motion. These parameters are similar for the two nebulae, pointing to similar physical conditions 
in these objects. The ro-vibrational cascade following Lya pumping is predicted to have low surface 
brightness signatures in the visible and near infrared. 

Subject headings: ISM:molecules — planetary nebulae individual (NGC 6853, NGC 3132) — ultra- 
violetilSM 



1. INTRODUCTION 



Molecular hydrogen emission is often seen in the in- 
frared (IR) observations of phot odissociation regions 
within planetary nebulae (PNe) ('Z uckerman fc GatlevI 
IT988: Kastner et al. 1996). Continuum pumped ul- 
traviolet (UV) fluorescence, as obse rved in reflection 
nebulae an d star forming regions l)Witt et al.l Il989t 
iFrance fc^ cCandliss 2005), has not been detected in 
PNe. The observed IR quadrupole emission origi- 
nates in elevated ro-vibrational levels of the ground 
electronic state The distribution of molecules 

in t hese levels is characteristic of either shock heat- 
ing ijZuckerman fc GatlevI 11988(1 or radiative cascade 
from excit ed electronic states l|Havashi et al.l 119851: 
iGatlev et~al . 1987). UV fluorescence can result in the 
dissociation of the mole cule 10 - 15% of the time 
llStecher fc WiniamsllT967ll . or a non-thermal population 
in the excited ro-vibrational levels of Shocks, 
with specific temperatures of ~ 2000 K, can thermally 
populate higher ro-vibrational states, re sulting in IR 
trans itions without the far-UV cascade l|Takami et al.l 
l200n|) . Detection of the UV fluorescence in conjunction 
with models of the infrared to ultraviolet scaling of H2 
emission can constrain the contribution of each mecha- 
nism to the IR flux. The line pumped H2 UV fluorescence 
described in this paper is consistent with measured IR 
line ratios and gas temperatur es that favor the shock 
excit ation scenario |Zu ckcrman fc GatlevI 1198^ IStorevI 
119841) . As shown in § 4.1, this type of UV emission, 
unlike continuum pumped fluorescence, has little effect 
on the observed IR spectrum. Weak specific features are 
expected in the optical. 

At typical diffuse ISM temperatures below 100 K 
(|Rachford et al.ll200^ . only H2 molecules in the v" = 
level of X^'E'^ make a significant contribution to the 
UV fluorescence. These molecules are excited by 
912 < A < 1110 A radiation into vibrational and ro- 
tational levels of higher electronic states (mainly -B^S^ 
and C^H„, Black & van Dishoeck 1987). As the temper- 
ature increases, higher ro-vibrational levels of A"^E^ be- 
come populated, and the threshold for fluorescent pump- 



ing moves to progressively larger wavelengths. Due to 
the long lifetime of the X^Y,'^ levels, the molecules can 
be further pumped to B^T,'^ or C^H„ by a sufficiently 
strong radiation field before they have time to radiate. 
This mechanism is predicted to be important for nebular 
gas densities in the range from 10^ < n < 10^ cm^'^ 
and a radiation field (in units of the interstellar average) 
of 103 < Go < 10-* feternbergl ll9891. 

We have found that in both NGC 6853 (M27, The 
DumbbeU Nebula) and NGC 3132 H2 UV fluorescence 
is dominated by lines pumped by nebular Lya from 
the excited v" = 2 level of X^T.+ to the 5^1]+ states. 

y 

Molecular hydrogen resonance fluorescence with Lya 
has been previously recognized as an important con- 
tributor to the UV flux in collisional environments 
of photodissoci ation regions. This effect was first 
pointed ou t by Shull (1978), and discussed further by 
■Black van Dishoeck (^19871 Lya pumping was pre- 
dicted to be significant mainly for T Tauri stars and 
Herbig-Haro objects, as has been subsequently observed 
(Brown et al. 1981; Schwartz 1983; Ravmond et al. 199'^ 
iHerczeg et al.l I2004f) . Other environments in which 
Lya pumping was d e tected include solar system objects 
llWolven e1^ all Il997t [Jordan et aIlll977D and accreting 
svstems l|Gizis et all 120051 iWood et al.ll2002D . In this 
paper we present for the first time evidence for resonant 
excitation by Lya of the B - X (1 - 2) P(5) 1216.07 A 
and S - X (1 - 2) R(6) 1215.73 A molecular hydrogen 
lines in planetray nebulae. The observed resonances and 
line ratios are a valuable diagnostic tool for molecular gas 
temperature, Lya line shape, and Doppler shift relative 
to the H2 absorber. 

The first object NGC 6853 is a 12,700 year old 
planetray ne bula l |(2^Dglj_^aljj2p02l) at a distance of 
417 ± 50 pc (iBen edict et al.ll2003l) . Its central star has 
a temperature 108,600 ± 6800 K l|Be'n cdict et al. 2003J. 
The distance to N GC 3132 is not well known, vary- 
ing from 1.6 3 kpc llTorres - Peimbert & PeimbertI ir977[l 
to 0.51 kpc l)PottaschrQ996ir while the temp erature of 
the i onizing star is estimated at 110,000 K l|Pottaschl 
^2^. The extinction, as measured by E{B — V) in 
this case, is of the order of 0.1 for the central star in 
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both cases (i Benedict et"al]l200llPottascb|[T99^ . How- 
ever, iMcCand liss et alJ l)200(il ^ iirgue that in the case of 
NGC 6853 the extinction is neghgible. In both cases, the 
detected Lya pumped H2 fluorescence requires a temper- 
ature of ~ 2000 K and column densities of ~ 10^* cm^^, 
as well as a Lya profile of ^ 0.4 A width, with a deep self- 
reversal. These parameters satisfy both short and long 
wavelength constraints, pointing to similar conditions of 
the molecular gas in the two nebulae. 

A description of the observations and data is found 
in Section 2. Data analysis follows in Section 3, and a 
discussion of the results in Section 4. 

2. OBSERVATIONS 

Four nebular observations of NGC 6853 were made by 
FUSE on 2004 May 26 using the low-resolution (LWRS) 
aperture (30"x30"; Figure Spectra were obtained in 
the 905 - 1187 A bandpass at a filled aperture resolution 
of ~ 0.33 A. A descrip t ion of the FUSE satellite can be 
found in IMoos et alJ l|2000() and on-orbit performance 
characteristics are described by ISahnow ct al. (2000). 
The data were acquired as part of the E120 guest ob- 
serving program. The pointings extend along the bright 
bar structure of the nebula, from the northeast to the 
southwest (Figure . Data for all four pointings were 
obtained in "time-tagged" (TTAG) mode, and processed 
using the CALFUSE pipeline, version 2.4.2. The average 
exposure time was 3374 seconds. The log of the observa- 
tions of NGC 6853 is presented in Table ^ 

We obtained optical spectra of NGC 6853 on 1999 
June 9 using the Double Imaging Spectrograph (DIS) 
on the Apache Point 3.5m telescope. The observations 
were taken with the low-resolution gratings, centered at 
8000 A and 4224 A with a dispersion of 7.0 A pixel"^ 
in the red side and 6.1 A pixel"^ in the blue side, and 
a resolution of roughly 2 pixels. The 300"x0'.'9 slit was 
centered at (19^ 59°" 36". 20, 22° 43' 01". 00), and ori- 
ented along the bright bar. Three slit positions (about 
2'.'7 total width) were averaged during the 900 s integra- 
tion. The observed mean Ha brightness along the slit 
was about 4.68 x 10"'^ ergs cm~^ s~^ sr~^, or 19416 R 
(conversion factor 632 Ahq(A) ). 

NGC 3132 was observed by HUT aboard Astro - 2 for 
946 seconds during orbital night on 1995 March 14. This 
spectrum, downloaded from the Multi-Mission Archive 
at the Space Telescope Science Institute (MAST), was 
acquired through the 10"x56" slit, enclosing a large re- 
gion of the 28" nebula, at an offset from the central stars. 
The exact slit coordinates are uncertain. The resolu- 
tion of the HUT spectrum is roughly 3 A. A descrip- 
tion of t he HUT instrumen t and data reduction can be 
found in iKruk et alJ l) 19951) . The signature of H2 emis- 
sion is present in the FUSE spectrum of NGC 3132 
(D1200401), but at a considerably lower signal-to-noise 
ratio. We chose to analyze the HUT spectrum because 
it allowes us to extend our study to longer wavelengths. 

3. DATA ANALYSIS 

3.1. NGC 6853 

The FUSE spectra of NGC 6853 - Position 3 are pre- 
sented in Figure 121 and Figure |3| The emission from the 
nebula is dominated by He II, C II, C III and molecu- 
lar hydrogen lines. The strongest lines observed in all 
four pointings are summarized in Tables 2-5. The lines 
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Fig. 1. — Optical image of NGC 6853 obtained during the 
commissioning of the FORSl instrument on the VLT 8.2-m. This 
three color image was made from a composite of B, [O III], and 
Hcf filters. FU SE LWRS aperture overlays are imposed showing 
the location of the four nebular pointings 1 to 4 from left to right. 
The optical image was obtained from the European Southern Ob- 
servataory. 



denoted by letters A — D are unidentified. Feature A 
is possibly nebular O III (962.425 A) at a slight veloc- 
ity offset. The strengths of the lines in the spectrum 
changes with the pointing. The highly ionized species 
dominate in the central regions and decline towards the 
outer shell. Molecular hydrogen lines are present in all 
four pointings, and are strongest in the third pointing, 
which coincides with the brightest optical feature. 

We explored the possibility that line-pumped H2 fluo- 
rescence is responsible for some of the unidentified lines 
we observe, computing models with a variety of nebular 
emission lines as the excitation source. Atomic emission 
hues of H I Ly-/3, Ly-(S, and Ly-7; He II A1215, A1085, 
A1025, A992, and A959; C II AA1334/35 and AA1036/37; 
the C III A1175 multiplet and C III A977; N I A1200; 
the N II A1085 multiplet; the N III A990 multiplet; and 
the O VI resonance doublet AA1032/38 were all con- 
sidered as possible H2 excitation mechanisms. None of 
the strongest predicted fluorescent emission lines are ob- 
served. 

The strongest H2 lines observed are excited by the 
B - X (1-2) P(5) transition at 1216.07 A in the red 
wing of Lya. Their wavelengths and branching ratios are 
listed in Tabled These transitions are sensitive to the 
Lya line width and Doppler shift. Although the closest 
resonance \s B — X (1-2) R(6) at 1215.73 A, Lyman- 
system fluorescence from v' = 1, J' = 7 is far weaker 
than expected. This can be achieved with a deep self- 
reversal of the Lya profile combined with a blueshifting 
of the molecular hydrogen lines with respect to the sys- 
temic Lya emission. The bueshift is necessary to restore 
to some extent the intensities of the lines pumped by the 
(1 - 2) R(6) resonance, which otherwise would fall in the 
center of the absorption trough (Figure 2)). 

NGC 6853, as well as NGC 3132, was modeled using 
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TABLE 1 
FUSE Observations of NGC 6853. 



Instrument Program 



Date 



RA 



Dec 



Exp. Time (s) 



FUSE - LWRS 
FUSE - LWRS 
FUSE - LWRS 
FUSE - LWRS 



E12001 
E12002 
E12003 
E12004 



2004-05-26 
2004-05-26 
2004-05-26 
2004-05-26 



19 59 44.96 
19 59 38.93 
19 59 32.35 
19 59 27.27 



-f22 44 50.1 
+22 44 04.0 
+22 42 35.1 
+22 42 10.3 



4176.0 
3071.0 
2778.0 
3470.0 




1130 1140 1150 1160 1170 1180 

o 




1130 1140 1150 1160 1170 1180 

Wavelength (A) 



Fig. 2. — FUSE spectra of NGC 6853 at the third pointing. Both the LiF lb (upper) and LiF 2a (lower) channels are shown to cover 
the whole 1115 - 1187 A interval. The H2 emission is prominent together with C II, C III and geocoronal He I. The Lya pumped H2 
fluorescence model is shown in red. 




1216.4 



Fig. 3. — The short wavelength FUSE SiC 2a spectrum of FiG. 4. — Exciting Lya profile absorbed by a H I column 

NGC 6853 at the third pointing. Strongest features are nebular density of 1 X lO^** cm"^ used for the NGC 6853 model. Velocity 

He II, C II, C III and geocoronal H I. offsets with respect to the molecular hydrogen reference system 

are 25 km s~^ for both emission and H I absorption. The H2 
absorbing transitions are indicated by dashed lines. The dotted 
line corresponds to an H2 absorption profile for a column density 
of 6 X lO^S cm-2. 



a fluorescent H2 emission code similar to the one de- 
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TABLE 2 
NGC 6853 Position 1. 



TABLE 4 
NGC 6853 Position 3. 



Line ID 


^obs 


FWHM 




Brightness 


Line ID 


^obs 


FWHM 


Brightness 




(A) 


(A) 


(10-6 


ergs cm"^ s"^ sr"^) 




(A) 


(A) 


(IQ— 6 ergs cm~^ 


C II 


946.01 


0.56 




1.83 ± 0.29 


He II 


933.40 


0.54 


1.21 ± 0.36 


He II 


958.52 


0.53 




0.35 ± 0.28 


He II 


937.30 


0.26 


0.91 ± 0.34 


N III + He II 


991.52 


0.21 




1.91 ± 0.39 


He II 


942.46 


0.30 


3.63 ± 0.44 


C II 


1036.88 


0.33 




1.37 ± 0.22 


C II 


946.03 


0.39 


4.75 ± 0.44 


N II + He II 


1085.50 


0.34 




25.75 ± 1.03 


He II 


949.30 


0.22 


2.06 ± 0.33 


C II 


1141.53 


0.26 




0.36 ± 0.13 


He II 


958.63 


0.29 


8.86 ± 0.54 


H2{1-1)R(3) 


1148.59 


0.43 




0.40 ± 0.14 


Un-ID (A) 


962.35 


0.27 


3.69 ± 0.45 


H2{1-1)P(5) + R(6) 


1161.76 


0.32 




0.71 ± 0.16 


N III + He II 


992.34 


0.32 


30.27 ± 1.02 


C III 


1175.53 


0.48 




1.97 ± 0.23 


S HI 


1015.42 


0.81 


1.20 ± 0.32 


H2(1-1)P(8) 


1183.06 


0.21 




0.92 ± 0.32 


S HI 
C II 


1021.16 
1036.81 


0.52 
0.33 


0.57 ± 0.20 
2.37 ± 0.29 


Note. — The C III line at 977.03 A 


is not listed due to contam- 


Un-ID (B) 
Un-ID (C) 


1062.34 
1065.79 


0.31 
0.30 


0.20 ± 0.25 
1.40 ± 0.26 


ination from scattered solar li| 


ght in the SiC channels. 


N II -1- He II 


1084.82 


0.34 


95.18 ± 2.18 












C II 


1141.48 


0.31 


1.50 ± 0.22 












H2(1-1)R(3) 


1148.48 


0.42 


0.68 ± 0.23 












H2(1-1)P(5) + R{6) 


1161.70 


0.42 


1.33 ± 0.22 




TABLE 3 






Un-ID (D) 


1167.44 


0.50 


1.65 ± 0.26 




NGC 68 


53 Position 2. 




C HI 


1175.48 


1.07 


28.33 ± 0.65 












. H2(1-1)P(8) 


1183.05 


0.26 


3.17 ± 0.49 



sr-i) 



Line ID 



X,t, FWHM 

(A) (A) 



Brightness 
-2 



(10 ergs cm s sr 



He II 


933.35 


0.26 


1.57 ± 0.43 


He II 


937.24 


0.21 


0.78 ± 0.37 


He II 


942.43 


0.29 


6.15 ± 0.50 


C II 


946.02 


0.46 


5.16 ± 0.48 


He II 


949.24 


0.43 


2.83 ± 0.38 


He II 


958.59 


0.29 


12.37 ± 0.61 


Un-ID (A) 


962.31 


0.27 


4.24 ± 0.48 


N III -1- He II 


992.30 


0.31 


48.83 ± 1.28 


S HI 


1015.48 


0.89 


1.65 ± 1.08 


S HI 


1021.32 


0.68 


0.72 ± 0.25 


C II 


1036.90 


0.42 


3.59 ± 0.34 


Un-ID (B) 


1062.62 


0.44 


0.31 ± 0.28 


Un-ID (C) 


1065.84 


0.47 


0.98 ± 0.27 


N II -1- He II 


1084.79 


0.32 


130.94 ± 2.46 


C II 


1141.51 


0.37 


1.26 ± 0.23 


H2(1-1)R(3) 


1148.63 


0.23 


0.62 ± 0.23 


H2(1-1)P(5) + R(6) 


1161.75 


0.40 


1.62 ± 0.23 


Un-ID (D) 


1167.30 


0.51 


1.24 ± 0.27 


C HI 


1175.53 


1.02 


45.55 ± 0.78 


H2(1-1)P(8) 


1183.26 


0.35 


2.01 ± 0.51 



TABLE 5 
NGC 6853 Position 4. 



Line ID 


^obs 


FWHM 


Brightness 




(A) 


(A) 


(10"" ergs cm~^ 


C II 


946.05 


0.38 


1.22 ± 0.30 


He II 


958.59 


0.24 


0.34 ± 0.32 


N III -1- He II 


991.63 


0.32 


1.56 ± 0.42 


C II 


1037.75 


0.42 


0.56 ± 0.22 


N II -1- He II 


1085.54 


0.34 


18.76 ± 0.97 


C II 


1141.52 


0.25 


0.21 ± 0.14 


H2(1-1)R(3) 


1148.50 


0.27 


0.12 ± 0.14 


H2(1-1)P(5) + R(6) 


1161.73 


0.21 


0.66 ± 0.17 


C HI 


1175.56 


0.29 


1.29 ± 0.23 


H2(1-1)P(8) 


1183.18 


0.31 


1.49 ± 0.34 



sr-i) 



scribed in iFrance et alJ l)200 5D. The exciting radiation 
field was assumed to be a Gaussian Lya emission line, 
with a FWHM of 0.4 A. In order to mimic self-reversal, 
absorption by neutral hydrogen with a column density 
of 1 X 10^'' cm"^ was used to modify the exciting Lya 
profile (Figure . As the nebular Lya line cannot be 
measured directly due to the opacity of the ISM and 
contamination from geocoroanal Lya, a total brightness 
was chosen as 2/3 of the Ha brightness of 19416.0 R 
measured using the DIS at APO (see § 2). According 
to the selection rules, 2/3 of the time the Ha transitions 
will result in an atom in the 2p state, which can then 
radiate a Lya photon during the 2p - Is transition. The 
other 1/3 of the time, the hydrogen atom will be left in 
the 2s state, unable to emit a Lya photon since a direct 
transition from this state to Is is forbidden in the dipole 
approximation (Spitzcr 1978). The total Lya brightness 
could be higher than derived from Ha measurements, 
due to H2 di ssociation mechanism s that produce H I in 
the 2p state l)Glass-Mauieaniri986j) . The continuum flux 



from the central star of NGC 6853 has not been included 
because continuum pumped emission lines are not de- 
tected in the nebular spectra. The likely explanation for 
the absence of continuum pumped fluorescence resides in 
the dominance of Lya photons over the 912 - 1110 A 
stellar continuum as well as the small covering factor in- 
volved in the absorption of continuum radiation, as will 
be discussed in § 4. 

The model parameters summarized in Table El were 
adjusted to reproduce the 1115 - 1187 A line strengths 
observed in the LiF IB channel. The LiF 2A chan- 
nel (Figure 0) shows a 20% decrease in the (1-1) R(3) 
brightness, which we attribute to calibration offsets. The 
molecular hydrogen temp erature of 2040 K was d erived 
from absorption spectra l)McCandliss et a L''2006). The 
H2 absorption is blueshifted with respect to the Lya pro- 
file by about -25 km s~^, value suppo rted by recent 
studies of NGC 6853 lIMcCandliss et al.ll200l . assum- 
ing Lya at the systemic velocity. The choosen total 
Lya brightness constrains the molecular hydrogen col- 
umn density of 6 x 10^^ cm~^. An H2 column density 
of few X 10^^ cm^^, as inferred from absorption spectra 
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TABLE 6 

Lyman system transitions pumped by Lyo out of the 

(v",J") LEVELS (2,5) AND (2,6). 



Line ID 


Wavelength (A) 


Branching Ratio (x 10*) 


(1— 0) R(3) 


lUyb. (O 


/D.bllu 


(1- 0) P(5) 


1109.31 


lUU.75z7 


(1— 1) R(3) 


1148.70 


osy.yoUb 


(1— 1) P(5) 


1161.82 


492.7111 


(1— 2) R(3) 




'7'70 QQO/1 

1 /5.oo24 


(1- 2) P(5) 


121D.U7 


926.0875 


(1- 3) R(3) 


1257. Bvi 


693.4022 


(1- 3) P(5) 


1271.90 


74z.zbbl 


(1- 4) R(3) 


1314.62 


1 /^o cr oner 

168.5805 


(1— 4) P(5) 


1329.14 


124.0794 


(1— 5) R(3) 


1 QVO A n 

lo/2.4y 


45. /uby 


(1- 5) P(5) 


loo7.o7 


123.9049 


(1— 6) R(3) 


1431.01 


ooO.o /o4 


(1- 6) P(5) 


1446.12 


8z5.4olz 


(1— ( ) tx[ o ) 


1 A CO V^l 

i4oy.o ( 


y44. / Uz4 


(1— 1) r[o) 


loU4. m 


ii/in /i/^c;c^ 

ii4y.4boo 


/I Q\ ■R?Q'l 

o J ixyo ) 




DD ( .OUDo 


( 1 q\ "D/" 
yl— o ) r\0 ) 




/ lo. / bzz 


(1— 9) R(3) 


IbUo.io 


01 '7 TOD y1 

21 / . /yo4 


[1— y ) r^[0 ) 


1 i^l 7 CO 

loi ( .oy 


iyb.b ( io 




1 ftc; c; OQ 

iboo.yvS 


OQ Qonn 
Zo.tSzUy 




1 ftRfl '71 


1 A CQflC^ 

i4.ooyo 


(1—11) R(3) 


1 rUo.oy 


0.0253 


l^i— ii ) r^[0 ) 


1 Vl fi 1 Q 


n nQ07 
U.UvJy / 


[L—IZ) ix[o) 


1 1 A A no 


U.Uoo4 


yi—LZ) r^yO) 


1 I^A i^l 
i ( 04. Oi 


n nyi OQ 
U.U4ZJ 


(1—13) R(3) 


1 HA A 1 

1 / /4.4/ 


n nnnn 

o.ooyo 


yL—Lo) t^yo) 


1 VQO 1 C 


n ni 1 n 
U.UiiU 


(1—14) R(3) 


1791.53 


n nnn/q 

O.UOOb 


(1— 0) R(6) 


iiuy.oD 


/?n noon 

by.y2ou 


(1— 0) P(8) 


1130.40 


00 A onn 

y5.4zoy 


(1— 1) R(6) 


1161.95 


Q'TT 00 ceo 


(1- 1) P(8) 


1183.31 


447.6342 


(1— 2) R(6) 


1215.73 


onn oono 

ouy.oo02 


(1— 2) P(8) 


1 OQ'7 QT 

Izo r .0 ( 


/ o.oOoz 


(1— 3) R(6) 


1 o'zi no 


/y / .ys / r 


(1— 3) P(8) 


1 onu 


Ton ^;nc;Q 
/2y. bOoo 


(1- 4) R(6) 


1o27.5d 


n /I T n/^ '7n 

247.0670 


(1- 4) P(8) 


1351.04 


"ion r* 

133.2556 


(1— 5) R(6) 


1385.01 


00 onno 

z2. ouyo 


(1- 5) P(8) 


1408.96 


111.5518 


(1— 6) R(6) 


1 A AO Q'7 


004.2o0( 


(1- 6) P(8) 


1467.08 


onn o/^'i^n 

800.3663 


(1- 7) R(6) 


1 50n 45 


1 009 0800 


(1- 7) P(8) 


1524.65 


1111.9502 


(1- 8) R(6) 


1556.87 


743.2822 


(1- 8) P(8) 


1580.67 


653.4937 


(1- 9) R(6) 


1610.95 


236.9584 


(1- 9) P(8) 


1633.84 


152.5215 


(1-10) R(6) 


1661.20 


20.8714 


(1-10) P(8) 


1682.49 


6.1068 


(1-11) R(6) 


1705.68 


0.0395 


(1-11) P(8) 


1724.43 


0.3495 


(1-12) R(6) 


1741.95 


0.0684 


(1-12) P(8) 


1756.90 


0.0004 


(1-13) R(6) 


1766.87 


0.0179 



ijMcCandliss et alJ l2006|) would require a much higher 
Lya brightness. 

3.2. NGC 3132 

In analyzing the NGC 3132 spectrum it is important 
to account for the background. The continuum has a 
very unusual shape, most likely due to the superposi- 
tion of the central stars and the inhomogeneous dust 
distribution within the nebula. A detailed analysis of 
the continuum is beyond the scope of this work, but a 



TABLE 7 
H2 Model Parameters. 



Parameter 


NGC 6853 


NGC 3132 


Excitation Source 


H I Lya 


H I Lya 


Lya Doppler shift 


25 km s-l 


30 km s-i 


Lya FWHM 


0.40 A 


0.45 A 


Lya total intensity 


12944 R 


217523 R 


T(H2) 


2040 K 


2040 K 


N(H2) 


6.0 X 10^* cm-2 


3.0 X lO^* cm-2 


N(H I) 


1.0 X lO^"* cm-2 


2.0 X 10i<^ cm-2 


b 


8 km s~i 


9 km s^^ 



TABLE 8 
NGC 3132. 



Line ID 


^obs 


FWHM 


Brightness 




(A) 


(A) 


(10~® ergs cm"'^ s^^ sr^^) 


H2(1-3)R(3) 


1259.19 


4.14 


8.62 ± 0.48 


H2(1-3)P(5) 


1273.38 


3.57 


16.81 ± 0.52 


H2(1-6)R(3) 


1431.73 


2.41 


11.86 ± 0.53 


H2(1-6)P(5) 


1446.67 


4.01 


26.75 ± 0.62 


H2(1-7)R(3) 


1489.85 


2.27 


23.12 ± 0.64 


H2(1-7)P(5) 


1505.84 


2.41 


22.59 ± 0.58 


H2(1-8)R(3) 


1548.16 


2.89 


15.54 ± 0.64 


H2(1-8)P(5) 


1563.10 


1.92 


14.35 ± 0.64 


H2(1-9)P(5) 


1617.69 


1.28 


7.04 ± 0.75 



rough baseline is needed for the model H2 spectrum. We 
set this background by manipulating the original HUT 
spectrum to extract a smooth continuum curve. A broad 
(width ~ 75 A) median filter was applied to the data at 
wavelengths between 1100 and 1750 A to remove most 
of the nebular and geocoronal emission lines. A boxcar 
smooth (width ~ 20 A) was then applied, leaving a back- 
ground continuum spectrum. Measured brightnesses for 
the strongest lines are given in Table|S| The observed H2 
emission lines provide a long wavelength confirmation of 
pumping by Lya photons. 

We use the NGC 6853 model as a starting point for 
NGC 3132 as a detailed study of the molecular struc- 
ture of NGC 3132 was not available. Previous obser- 
vations and models of NGC 3132 (Montciro ct al. 200(| 
iBassgen et al.l 119901) estimate the nebular Ha luminos- 
ity at few X 10'*'' photons s~^, too low to predict 
the observed level of H2 fluorescence. This value was 
derived using a distance of 670 pc, which might be an 
underestimate. The estimated total Ha luminosity of 
NGC 6853, approximated as an elliptical nebula with 
semimajor axes of 2.'5 and 4', at a distance of 417 pc, 
is 8.55 X 10''^ photons s~^ Using for the NGC 3132 
0.'45x0.'7 nebula a total Ha luminosity of 1 x 10'*'' pho- 
tons s~^, closer to the value for NGC 6853, we obtain a 
reasonable fit with column densities of 3 x 10** cm~^ 
for the molecular component and 2 x 10*^ cm~^ for 
the atomic hydrogen. The molecular hydrogen absorp- 
tion is blueshifted with respect to the nebular Lya by 
-30 km s~*, in agre ement with other measurements 
l|Monteiro et all 12000(1 . A summary of the parameters 
used is given in Table [7| The overprediction of line 
strengths in the 1200 < A < 1300 A region (Figure [H 
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and shortward of 1200 A (not shown) is caused by an in- 
complete removal of the combination of geocoronal Lya 
and the (1-3) H2 emission band together with opacity 
effects. These shorter wavelength lines connect to more 
populated lowe r vibrational levels (Table Et, leading to 
self-absorption l|Herczeg et al.ll2004t IWood et al.ll2002(l . 

4. DISCUSSION 

Although absorpti on was observed in both nebulae 
ijMcCandliss et alJ2006tlSterhng et al.l2002() . continuum 
pumped molecular hydrogen emission was not detected. 
Predicted continuum pumped UV fluorescence for an H2 
column density of 7.9 x 10^^ cm~^, me asured from ab- 
sorption spectra ijMcCandhss et al.l2006(l . is estimated at 
a level of ~0.7 x 10~^ ergs cm~^ s~^ sr~^ A~^, consis- 
tent with the noise level in Figure[21 However, the larger 
column density used in the Lya pumping model, would 
predict a much higher continuum fluorescence, which 
is not observed. One possible explanation comes from 
the dominance of Lyman continuum photons over the 
912 - 1110 A photons which contribute to continuum 
H2 pumping. In the central H II region of a PN, Ly- 
man continuum photons are readily converted into Lya 
line photons, following a series of ionizations, recombi- 
nations and cascades. Integrating a stellar model with 
T,ff = 110 OOP K, log(5) = 6.7, and = 0.56 Mq 
l|Rauchll2003t) . we find that there are over 16 times more 
photons cm~^ emitted in the Lyman continuum (5 - 
911.7 A) than in the 912 - 1110 A region. Additionally, 
the Lyman continuum photons are concentrated into an 
emission line (FWHMli,q = 0.4 A), compared with the 
far-UV stellar continuum that spans roughly 200 A. For 
an H II region with an electron temperature of 12,000 K 
(i.e.- NGC 6853; Pottasch et al. (1982)) the recombina- 
tion efficiency for generating Lya photons is roughly 67% 
()Spitzeilll978|) . which gives us 5.5 x 10^^ photons s ^. 
The total Lya luminosity obtained this way is still low 



compared to the value of 8.55 x 10^^ photons s~^ de- 
rived from the Ha brightness. This shows that the ratio 
of Lya photons to 912 - 1110 A continuum photons might 
be even higher than estimated. 

However, the dominance of Lya photons over the 912 - 
1110 A stellar continuum alone is not sufficient to com- 
pensate for the larger ground state populations and 
higher oscillator strengths contributing to other Lyman 
and Werner transitions out of v" = 0. Assuming that 
molecules reside in high density globules, the covering 
factor involved in continuum excitation is much smaller 
than for the diffuse radiation field of scattered Lya. If we 
let Bo be the surface brightness of the exciting field, the 
total flux absorbed by an H2 clump will be Bpfi, where fl 
is the solid angle subtended by the absorber, as seen from 
the source. The redistributed brightness is then radiated 
into 47r, so that we can define an effective surface bright- 
ness seen by the absorber as (Boil) / 47r, where ^ / 4it 
is the covering factor. For a 10" x 10" globule at a 50" 
separation from the central star, the continuum covering 
factor is about 0.0032, while for the nebular Lya it is 
likely to be unity. This estimate takes into account that, 
while the continuum photons are coming mainly from the 
star, the Lya photons are produced and scattered in the 
nebula, so that the globules are effectively embedded in 
a diffuse Lya radiation field covering all Att sr. 

The broad (98.6 km s"^ FWHM) Lya line required to 
pump the H2 lines in NGC 6853 is likely to originate in 
the complex shell structure of the nebula. As shown by 
Mcaburn et al. (2005), in the case of Ha the wide profile 
is a result of combining the motions of the outer shell at 
35 km s""'^, the inner shell at 13 km s~^ and the bulk 
motion less than 7 km s^^ in the central He II region. 
We infer a similarly broad velocity structure should be 
present in NGC 3132 to produce a Lya of comparable 
width. A narrower velocity distribution would be al- 
lowed if we take into account the redistribution of the 
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Lya photons in the hne wings for large optical depths. 

IR observations of NGC 6853 and NGC 313 2 made to 
date l)Storevlll98^ iZuckerman fc Gatlevlll988l) conclude 
that the H2 spectrum is shock excited based on the ratios 
of the S(l) (1-0) and S(l) (2-1) lines. Moreover, this con- 
clusion is supported by the non-detection of continuum 
pumped UV fluorescence, both in t hese spectra and i n 
rocket observations of th e Dumbbell llMcCandhssll200lD 
Howe ver, recent studies lIHora et al]ll999t iTakami et al.l 
12000(1 show that in dense enough environments the ratio 
of the S(l) (1-0) and S(l) (2-1) hnes can appear thermal, 
even if the IR emission is excited via UV pumping. In 
these cases, the study of transitions from higher energy 
levels becomes important in order to distinguish between 
the two mechanisms. The detection of II2 continuum ab- 
sorption towards the central star of NGC 6853 shows that 
continuum pumped fluorescence takes place, although at 
a level allowed by a non-detection in the FUSE spec- 
tra. This by itself does not rule out a partial UV contin- 
uum pumping of the IR lines. However, the thermal pro- 
cesses are thought to be dominant, since absorption out 
of vibrational levels v" > 2 is not observed, in contrast 
with reflec tion nebulae where fourescence is important 
l)Mever et al. 2001). The possibility of shock excitation 
in NGC 6853 and NGC 3132 is also revealed here by 
the high temperature required by the prese nce of Lyct 
pumping. More in-dep th analysis of PNe l(Hora et al.l 
119991 iDavis et alJl200 reveal that in most cases where 
both the rotational and vibrational temperature of H2 
are around 2000 K, shock heating is the dominant ex- 
citation mechanism of IR lines. In order to confirm the 
shock scenario, we need a better understanding of the gas 
motions in the two objects, in addition to observations 
shortward of 2 /im. 

The observed presence of both ionized species and 
molecules in the same aperture, and the small cov- 
ering factor necessary to explain the non-detection of 
continuum pumped lines confirms the current view 
of molecules residing in globules surrounded by an 
ionized medium, the source of the exciting Lya ra- 
diation llMc aburn & Lopez 1993: Huggins et al. 2002; 
ISpeck et ai.H2003ll . The presence of globules helps the 
survival of hydrogen molecules, shielded from photodis- 
sociation by hard stellar radiation. Overall, the H2 col- 
umn density derived for the Ly alpha pumping model is 
much higher than the one derived from absorption spec- 
tra (10^* cm~^ vs 10^^ cm~^), suggesting different prop- 
erties of the regions they are probing. Taking into ac- 
count that the globule size in NGC 6853 is likely to be 
less than 10" , the FUSE observations average between 
high and low density areas, where the low density re- 
gions are thought to have a negligible extinction. The 
low slit filling fraction for the Lya pumped radiation in 
both FU SE and HUT observations allows us to reconcile 
the derived value of fewxlO^* cm~^ for the H2 column 
density with coulmns of 10^^ cm^^ sugges ted for dense 
globules fr om visual exctinctio n arguments ijBohlin et al.l 
Il978: Meab^irn fc T.one7ilT99l . 

4.1. Ro-vibrational Cascade Modified by Lya Pumping 

We constructed the ro-vibrational spectrum as a com- 
bination of a non-equilibrium cascade following Lya 
pumping and a steady state emission due to the equilib- 
rium distribution of molecules on the X^Y.'^ state levels. 



For the non-equilibrium cascade, we start with the rates 
at which the X^'E'^ levels are populated following Lya 
pumping, and evaluate the line strengths according to the 
branching ratios. CoUisional effects are ignored for the 
purpose of this estimate. As a proxy for the equilibrium 
populations for the steady state case we use the ground 
state populations measured from absorption spectra and 
extrapolated to higher vibrational and rotational levels 
using a 2040 K thermal distribution with a total H2 col- 
umn density of 7.9 x 10^^ cm~^. For this reason here- 
after we will refer to this steady sta te model as ther- 
mal. T he transition probabilities from iWolniewicz et all 
l)1998D are used to derive the output spectrum, shown 
in Figure El The solid black line contains contributions 
from both thermal and fluorescent pumping, convolved 
with a 15 A Gaussian. Overplotted in red is the thermal 
contribution alone. The lines showing the largest con- 
tribution from fluorescent pumping are numbered, and 
their principal components are listed in Table [U] The 
line strengths are given as upper limits and will be af- 
fected by coUisional de-excitation. CoUisional excitation 
is not important for these high lying vibrational states. 
Using for the thermal model the same total H2 column 
density of 6 x 10"'^^ cm~^ used for the Lya pumping, the 
lines with the highest contribution from fluorescent cas- 
cade wiU be 1,2,3,4,8,9 and 10, as numbered in FigureEJ 

The ro-vibrational cascade shows most of the Lya 
pumping specific features in the visible and near-IR part 
of the spectrum. We find little deviation from a pure 
thermal emission longward of ~ 1 /im. As a consequence, 
line pumped UV fluorescence might be present even when 
the measured IR line ratios around 2 /im are consistent 
with a thermal distribution of the H2 molecules in the 
ground electronic state. 

The detection of these lines represents an observational 
challenge from the ground. In Figure [S] the positions of 
the strongest predicted lines are indicated on the air- 
glow emiss ion spectru m reconstructed from UVES ob- 
servations l)Hanuschikl 12003) at a resolution of 15 A to 
match the resolution of our model. The lack of air- 
glow features in the 0.577 - 0.583 fim interval is due 
to the chip gap in th e UVES spectrum. However, 
lOsterbrock et al] l|1996|) shows that there are no lines 
identified in this region. The ro-vibrational cascade at a 
level of few x 10~^ ergs cm~^ s~^ sr"-'^ fj,m~^ is about 
1000 times smaller than the typical airglow lines and 
contiuum in the 0.56 - 1.015 /zm interval. A much 
lower background level and a higher spectral resolution 
(~ few A) may allow the detection of the (9-3) S(4) line 
at 5810.36 A and the (8-3) S(6) line at 6681.40 A (Fig- 
ure EI). 

Lya pumping redistributes molecules from the (v",J") 
levels (2,6) and (2,5) among the ro-vibrational levels of 
the ground state. However, at the level of the observed 
UV fluorescence, a significant deviation from the ther- 
mal populations of the (v",J") levels (2,6) and (2,5) 
is not found. Correcting for the ro-vibrational tran- 
sitions that repopulate these levels, we find that the 
population decrease relative to the thermal distribution 
(A Nyj/Nlf™"^) is about 0.0107 and 0.0012 for the 
(2,6) and (2,5) levels respectively. Wc find more signifi- 
cant deviations from a thermal population among the ro- 
tational levels of the v"=0 state. These are likely to be af- 
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Fig. 6. — Predicted ro-vibrational cascade following Lyo pumping in NGC 6853 (see text). The red dashed line is the ro-vibrational 
cascade model for a thermal population at 2040 K. The black continuous line has both contributions added. The main contributors to the 
numbered lines are listed in Tabled 
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Fig. 7. — Airglow spectrum reconstructed from UVES data 
IHanuschikl l2003) . with a resolution of 15 A. The positions of the 
strongest lines following Lyg pumping are indicated by the green 
lines, numbered as in Figure 161 



NGC 6853 and NCG 3132. The observed Ha fluores- 
cence gives us valuable information about the conditions 
of the radiation field and dynamics of the molecular gas 
in planetary nebulae. The input parameters are similar 
for both objects, suggesting that we do not see an iso- 
lated phenomenon. The UV H2 spectrum is likely due 
to thermally excited molecular hydrogen, shielded from 
the UV continuum inside globules surrounded by strong 
nebular Lya emission. We emphasize that while molec- 
ular hydrogen pumping by Lya shows strong features in 
UV spectra, it could remain unobservable in the IR and 
visible. UV observations are thus a requirement for in- 
vestigating this process in various environments where 
hot H2 is exposed to the presence of Lya. 



fected also by collisional redistribution and do not match 
the deviations measured from UV absorption spectra 
(j^cCandliss et al. 2006). Deviations in the populations 
of higher vibrational levels (v" > 2) are predicted to re- 
sult in a specific signature in the ro-vibrational spectrum 
(Figure El. 

5. CONCLUSION 

Line pumping by Lya is required to qualitatively ex- 
plain the H2 emission features in the far UV spectra of 
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RO-VIBRATIONAL CASCADE, THERMAL CONTRIBUTION NOT 
INCLUDED. 
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